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51 Derivation of MSD for non-motile bacteria 

Integrating equation (3) gives 

x(0 = C ! f €(u)du- (SI) 
Jo 

The mean squared displacement (MSD) is defined by 

(||x(0H 2 ) = C 2 ( f f t(u)-t(v)dudv) = 2{; 2 f f (Uu)Uv)) dudv. (S2) 
\Jo Jo / Jo Jo 

The second equality is obtained by noting that the components g x (i), g y (t) of the 2D vector f (t) 
are independent and identically distributed and applying Fubini's theorem [? ]. Substituting (2) 
and integrating, we obtain 

kT 

(\\x(t)\\ 2 ) = A—t = 4D t t, (S3) 

it 

where D t = kT/£, is the translational diffusion coefficient. Equation (S3) is identical to results 
presented in the literature on this subject [? ]. 

52 Modified censoring method 

Figure S 1 shows the analogous plots to those in Figure 8 of the main text, computed using the 
non-chemotactic dataset with the same censoring parameters as used in the original study [? ]. 
The departure from linearity in the MSAC is more significant at shorter timescales than those 
apparent with the modified censoring process. Even more striking is the departure from normal- 
ity in the observed distributions of framewise angle changes. This discrepancy is significantly 
reduced when the modified censoring process is used, suggesting that the 85 additional tracks 
removed from the non-chemotactic dataset in our modified censoring process are from a dis- 
tinct subpopulation with greatly altered motility. This is further supported by the appearance of 
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the additional tracks, ten of which are plotted in Figure S2. Comparing Figure S2 with Figure 
1 in the main text, we see that these tracks are more tortuous than the majority of tracks in 
this dataset, which appear smooth and contain no stops. It is likely that these tracks arise from 
bacteria with damaged motility apparatus, or collections of two or more attached bacteria. 



S3 Derivation of observed cell displacement distribution 

Starting from equation (19), we may write 

x(0 = C 1 f ?(«)d«=^B(0, (S4) 
Jo 

where B(t) = (B x (t), B y {t)) denotes the 2D Wiener process [? ], whose components are indepen- 
dent and normally distributed with zero mean and variance t. The multiplicative factors on the 
right-hand side of (S4) are required to correct for the non-unit variance of the white noise force 
%(f). The characteristic functions of B x {t), B y {t) are given by 

= ( e *By(f)) = e -" 2 " 2 _ (S5) 

The absolute displacement R(t) over a time step r = t - f for t > t' is given by the square root 
of a sum of squared normal random variables, 



R(t) = \\x(t) - x(t')\\ = ^ X { t )-x{f)f + {y{t)-y{t')) 2 . (S6) 

Since x(t) is normally distributed, so too is the difference, x(t) - x(f), and thus R(t) has proba- 
bility density function (PDF) 



— f(— 
V2Ar \ V2Ar' 



f R (r; r) = - 1 ==f x I - 1 ==, 2 1 , (S7) 



where f x (-, k) denotes the PDF of the^ distribution with k degrees of freedom. 
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S4 Derivation of analytic results for run-only model 

We begin our analysis of the run-only model (8)-(10) by deriving expressions for the first and 
second moments of the angle change, <p{t). Integrating equation (10) subject to the initial con- 
dition 0(0) = 0 and using (S4), we obtain 

0(0 = Tj2D r B x (t\ (S8) 

where D r = kT/£ r denotes the rotational diffusion coefficient. This, rather than £ r and T, is the 
quantity of interest in most experimental studies in this area. Using properties of the Wiener 
process, we see that the mean angle change is simply {(f>{t)) = 0, while the mean squared angle 
change (MSAC) is given by 

<0 2 (O> = 2D r t. (S9) 

We now compute the first and second moments of the position vector. Substituting (S8) into 
(8), integrating subject to x(0) = 0 and rewriting the cosine function in complex form yields 

^i^2D r B(s) +e -i^mB(s^ ds (S1()) 

Taking the mean of both sides, applying Fubini's theorem [? ], substituting (S5) and integrating, 
we obtain 

(x(f)) = £-(l-e- D ' t ). (Sll) 

From (Sll) we see that (x(t)) « ct for small t, indicating that the propulsion of the bacterium 
initially dominates over rotational diffusion and the motion is approximately ballistic. Over 
a long timescale, (x(t)) tends to the constant ratio c/D r , indicating that rotational diffusion 
overrides the effect of the propulsion and that, on average, our model bacterium can only travel 
a certain distance from its starting point before being constrained by the randomising rotational 
buffeting. Since we assume that y(0) = 0(0) = 0, by symmetry we have (y(t)) = 0 for all t > 0. 
The second moment of x(t) is given by 

(x 2 (t)) = c 2 ^ cos ( yj2D r B(u)) du ^ cos ( V2£>>(v)) dvj . (S 1 2) 



x(t) = c \ 
Jo 
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Interchanging the order of integration using Fubini's Theorem and using standard properties of 
the Wiener process we obtain 

A similar process may be followed to obtain the second moment in the y-coordinate, which is 
given by 

^-tefS^-i^-iH- <S14) 

From (S 1 3) and (S 14), we see that on a timescale that is short relative to the diffusion coefficient, 
t «; l/D r , we have (x 2 (t)) ~ c 2 t 2 and (y 2 (t)) « c 2 t 2 /2, which again indicates that movement is 
ballistic. This is also the case for the VJ process over short timescales [? ]. For longer times 
t » l/D r , both second moments vary approximately as c 2 tlD r . This asymptotically linear 
behaviour of the second moment is characteristic of a diffusive process and is again similar to 
that observed in the VJ process. Coupling this result with equation (SI 1), we see that the mean 
position of the bacterium reaches an asymptotic value, but the bacterium continues to explore 
its environment in a diffusive manner. Furthermore, we may combine these equations to obtain 
an expression for the standard deviations of x(t) and y(t) using the standard result 

std(x(0) = VV(0> - (x(t)) 2 . (S15) 

We obtain 

std(x(0) = ^ (^e- 4D >< - e- 2 ^ + + D r t - ^ , (S16) 

std(y(0) = (-T^- 4ZV + \e~ Drt + D r t - ^ . (S17) 



S5 Derivation of analytic results for run-and-stop model 

We begin our analysis of the run-and-stop model (1 1)— (13) by noting that the solution to equa- 
tion (13) matches that of (10). We now consider the stochastic switching force F(t) in equations 
(1 1)— (12). Let p(t; s) = P(F(t) = l\F(s) = 1) denote the probability that the bacterium is run- 
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ning at time t, given that it is running at time s < t. The evolution of pit; s) is governed by the 
forward equation [? ] 

^=A r -(A r + A s )p (S18) 

subject to pis; s) = 1, hence we have 

p(t; s) = — -(/t r + A s e-^*- s) ). (S19) 

A S + Ay 

Since F = 1 is the only allowable non-zero value, the ensemble average of F(t) is given by 

1 



(F(t)) = p(t; 0) = — - -iA r + A s e^ M ) (S20) 

A s + A r 

and its autocorrelation function is given by [? ] 

(F(s)F(t)) = (F(t)\F(s) = 1><F(j)> (S21) 
= pit;s)pis;Q) (S22) 

= (A r + A s e- (As+Ar){ '- s) ) (A r + A s e~ (AM ) , (S23) 

\A S + Ay J 

for 0 < s < t. Integrating (11) gives 

xit) = c ^ Fis) cos ( ^JlDyBis)) ds. (S24) 

Noting the independence of Fit) and Bit), we carry out similar manipulations to those used in 
deriving (Sll) to show that the ensemble average of xit) is given by 



(xit)) = ' 



Ag + Ay 



A(l _ e -D rt) + (1 _ e -(D,- + A s+ A r)t) 

Dy Dy + Ag + Ay 



(S25) 



As for the run-only model, we have (yit)) = 0 by symmetry. The second moment of xit) is 
given by 

(x 2 it)) = c 2 J J (F(m)F(v)) ((cos [ V2A^(5(v) + £<»)]) + (cos [ V^CW - £<»)])) du dv. 

(S26) 
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Substituting equation (S23) into (S26) and using standard properties of the Wiener process gives 

(S27) 

The integral in equation (S27) is straightforward to calculate but is not presented here for the 
sake of brevity. A similar expression may be derived for the second moment of y(t), and is given 
by 

(/(,)) = |_^!_J 2 ^ ^{e' DAv ~ u) - e -^ (3u+v) J|l + ^ e - (i ^ )H J|l + ^e- (AMv - u) y u <\v. 

(S28) 

Expressions for the standard deviations of x(t) and y(t) may be then be computed from (S27) 
and (S28). 

In the case of the run-and-active-stop model, the moments can be derived using similar 
steps to those shown above. However, we do not solve for the moments of this model here as 
the result would be algebraically cumbersome and not instructive for our purposes. 



Ay 



i , 4? -a+^xv-w) 



du dv. 
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Figure SI. Extracting characteristics of noise from the R. sphaeroides non-chemotactic dataset 
obtained by ? ] with no modification to the censoring process, (a) The observed MS AC (black 
crosses) with linear fit (dashed line). The slope of the fit is 0.255 rad 2 s _1 . (b) The estimated 
PDF of angle changes (solid line), overlaid with the wrapped normal distribution with variance 
computed from the data (dashed line). The three shades denote different values for the 
sampling interval, r. Black: r = 0.02 s; dark grey: r = 0.1 s; light grey: r = 0.2 s. 
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Figure S2. Ten randomly-selected tracks from the non-chemotactic dataset that are included in 
the original study [? ] but excluded with the modified censoring process used here. 
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Figure S3, (a) 20 sample trajectories generating by simulating the run-only model (8)— (10). 
(b) Mean (solid line) and mean + standard deviation (dashed lines) for x (black) and y (grey) 
position of a particle in the run-only model. 
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Figure S4. (a) 20 sample trajectories computed by solving the run-and-stop model (11 )— (13). 
(b) Mean (solid line) and mean + standard deviation (dashed lines) for x (black) and y (grey) 
coordinates of a particle in the run-and-stop model. 
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Figure S5. Histograms of the x and y position of simulated particles in the run-only model 
(8)— (10) with sampling time step r = 1 s. 



